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Abstract. The cellular automata discrete dynamical system is consid- 
ered as the two-stage process: the majority rule for the change in the 
automata state and the rule for the change in topological relations be- 
tween automata. The influence of changing topology to the cooperative 
phenomena, namely zero-temperature ferromagnetic phase transition, is 
observed. 



1 Introduction 

Lattice models are the basic instrument in the study of phase transitions in 
equilibrium statistical mechanics (see, e.g., p^). The Ising model is the sim- 
plest model of nearest-neighbor ferromagnetic interactions where the collective 
features are studied. Traditionally, the phase transition has been considered in 
regular lattices. According to the renormalization group theory, detailed struc- 
ture of interactions as well as the structure of network connections are irrelevant. 
Therefore, for example, cellular automata have been continuously tested with the 
hope that these systems can imitate nature [2]. 

If the cellular automata are considered on a regular lattice and each site i 
of the lattice is occupied by a binary variable, s,, which assumes the values ±1, 
then referring to these variables as spins which point either 'up' (+1) or 'down' 
(—1), we employ the terminology of the Ising model. The Ising nearest-neighbor 
interaction in case of cellular automata becomes usually a simple majority rule, 
it is: if most of the nearest neighbors of i spin point up (down) at time t then 
the value of the i spin at time t + 1 is up (down). This deterministic rule is 
interpreted as the zero-temperature approximation to the Ising interactions. The 
research here is aimed on the ergodicity, i.e., on the uniqueness of the final state. 
Especially, one asks if starting from a random distribution of spins states one 
ends up with a homogeneous states of all spins up or all spins down. To mimic the 
temperature effects one adds a stochastic noise to perturb with some probability 
the execution of the deterministic dynamics. The research here is focused on 
properties of the transition from the non-ordered phase to the ordered phase - 
ferromagnetic phase transition. 

Recent research in the structure and topology of real networks [.'il lj has shown 
that social, biological, technological networks are far from being regular. How- 
ever, they also are being far from a random network. It has been shown that 



networks called small-world networks and scale-free networks exhibit mixed prop- 
erties of regular lattices and a random graph |3I5| . This has trigged to study 
of standard models of statistical mechanics in these networks. The transition 
from non-ordered phase to the ordered phase in the spin system spanned on 
the small-world network shows a change in behavior from the regular case to 
the mean- field characteristics |7I12| . Interesting that, depending on the network 
construction, the phase transition at the finite temperature exists or not UJ. 
Moreover, when the scale-free network is employed then the critical properties 
of the Ising model are different from those ones observed in the regular networks 
|13I14| . Many authors have considered other problems in these new networks: 
percolation properties jS], the spread of infection diseases [§] social relations |10j . 
computer organization [TT] . 

Small-world networks are intermediates between regular lattices and the ran- 
dom graph. A small-world network is generated by rewiring with a probability p 
the links of a regular lattice by long-distance random links [5]. The presence of 
a small fraction of 'short-cuts' connecting otherwise distant points, drastically 
reduces the average shortest distance between any pair of nodes in network. 

In case of cellular automata such the small-world network means that the 
cellular automata is no longer homogeneous — inhomogeneity is due to the 
topology. Sets of nearest-neighbors differ from a site to a site. 

There is one important property which is ubiquitous in biological and artifi- 
cial networks and which is missed in the small- world network — the distribution 
of the vertex degree does not show wings decaying as a power-law. Two important 
ingredients have been shown to be sufficient to generate such feature: growing 
number of vertices and preferential attachment of links. The well established 
Barabasi- Albert model based on these two mechanisms has been proposed [3]. 
The distributions of the vertex degree in this network is of the form P(k) oc fc~ 7 
where 7 6 [2,3] and k a vertex degree. Thus adapting the Barabasi- Albert idea 
to cellular automata one obtains a computing system with the topology that 
evolves. 

In this paper we address the question of the role played by the topology 
in the zero-temperature ferromagnetic transition considered in the cellular au- 
tomata of spins. The starting system is the cellular automata on a square lattice. 
The square lattice used here is strongly clustered in the sense that if j and k 
are neighbors of i then there is a short path between them that does not pass 
through i. In order to keep the high clustered property (typical for the small- 
world network) and reproduce the distribution of the vertex degree with the 
power-low decay (sign of the scale-free network) we propose the evolution rule 
for the network. Hence, the cellular automata time step means the change of the 
spin states and the change of the network topology. It will be shown that with 
increasing p- the rewiring parameter, the cellular automata work as the solver 
to the density classification task |15j . highly accurately, namely our cellular au- 
tomata converge to a fixed point of all l's if the initial configuration contains 
more l's than -l's by 2%, and, symmetrically, to a fixed point of all -l's if the 
initial configuration contains more -l's than l's by 2%. 



2 The model description 



We present the model of the spin cellular automata located in N vertices among 
which there is an explicit integer parameter Nq fixing the total number of edges. 
These N spins are located to form a square lattice with the periodic boundary 
conditions. Hence, if the linear size of the lattice is L then N = L x L and 
No = 4iV. We start with the ordinary nearest-neighbor relations, i.e., for any i 
site the nearest-neighbor spins are located in the following set of vertex indices: 
N(i, 0) = {i — L, i — 1, i, % + 1, i + L}. Initially, each spin state is set up randomly 
with the probability p. 

The evolution step consists of the two subsequent steps: (A) the asynchronous 
stochastic evolution of topological relations and (B) the synchronous determin- 
istic majority rule. 



2.1 Evolution of the lattice 

In each time-step t and for each vertex i the following rearrangement of the edges 
is performed: from the set of neighbors N(i, t) of the vertex i (at the given time 
t) one element j is chosen at random with the probability p. Thus p relates to 
the parameter of the small- world evolution . Then, another vertex k is picked 
up from the set of all vertices N — {i} randomly. A new edge is created between i 
and k independently of the fact that there was or not an edge between them. The 
edge between i and j is deleted. Hence the total number of edges Nq is conserved. 
This evolution we will call stochastic on the contrary to the evolution in which 
some preferences are applied. The proposed preferences are in the agreement 
with the basic conviction how to obtain a scale free distribution for the vertex 
degree: 'the richer you are, the richer you get' — principle 

The preferences are as follows: in each time step t and for each vertex i, from 
the set of neighbors N(i,t) an element j is chosen with the probability modified 
by the degree of j vertex. Namely, the probability that j neighbor is unlinked is 

PI 1 — ] I 1 ) 

where deg(j) denotes the degree of j vertex, T denotes the value of the unlink 
threshold. The degree of a vertex denotes the number of edges attached to the 
vertex. By Q if deg(j) > T then it is impossible to unlink j. This process will 
be called the intensional detachment. The edge between i and j is deleted. Then 
a new vertex is randomly chosen but with the preference to link to a vertex with 
the high vertex degree. It is, a randomly picked up vertex k is accepted with 
probability 

deg(k) ^ 

Thus if deg(k) > T then randomly chosen k is certainly accepted. This process 
will be called the preferential attachment. 

In the following we will name: 
model : the completely stochastic edge evolution, 



model 1 : the intcnsional detachment and the stochastic attachment, 
model 2 : the stochastic detachment and the preferential attachment, 
model 3 : the intensional detachment and the preferential attachment. 

2.2 Evolution of spins 

Synchronously, in each time step, the majority rule is applied to set the state of 
every spin in the next time-step. The rule acts on all spins belonging to the the 
set of nearest neighbors N(i, t) of the given ith spin. In the case when the result 
of voting is zero, the ith spin does not change the state. 

3 Results 

The system simulated is a square lattice with L = 200 and periodic boundary 
conditions. We start with a random distribution of spin states with probability 
p to set the spin state to +1. In each time step we measure magnetization - the 
sum of all spins states normalized by L 2 . If the magnetization is 1 or the time 
step reaches the limit of the hundredth step, then we stop the evolution and 
record the following characteristics: magnetization, susceptibility (the variation 
of the magnetization), vertex degree and number of time steps to reach the 
stable state. The experiment is repeated 200 times for each p and p value. The 
threshold value for the preferences, see JIJ, © , is T = 8 in all experiments. 

3.1 Final state characteristic 

The sharp change in properties of the final state of cellular automata is ob- 
served when p is crossing 0.5, see FigQ] Fig[3 The final state switches between 
the state of all spins down if p < 0.5 (magnetization = —1) to the state of 
all spins up if p > 0.5 (magnetization — 1), Fig^ Depending on the model of 
the network evolution this transformation goes fast: model 0, or less rapidly: 
model 3. The transformation can be estimated by measuring the intervals of p 
around the 0.5- value which lead to the final states with with magnetization 
close to ±1, see Fig.l. More accurate marks for these intervals can be found by 
observing susceptibility, see Fig. 2. In case of unperturbed network (p = 0) the 
susceptibility takes the parabola-like shape on the log-plots. With increasing p 
the suscpetibility becomes the sharply pointed curve with maximum at p = 0.5. 
The basic step in p in the experiments performed is Ap = 0.02. Therefore the 
peaks of susceptibilities are recognized within the interval (0.48,0.52). In the 
small windows of Fig. 2 one can observe how fast and which way the parabola 
transforms into the sharply pointed curves. The significant qualitative change 
occurs when p changes from p = 0.01 to p = 0.02. In particular, if one assumes 
that the interval of critical changes means the susceptibility greater than 1 (at 
average this condition is equivalent to the demand of the absolute magnetization 
being smaller than 0.9), then one finds: 

— p S (0.30,0.70) in case of cellular automata without the network evolution 
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Fig. 1. Magnetization of final configurations vs initial probability for a spin to 
be UP (p) for different values of short-cut edges (p) what is indicated by curve's 
labels: (A) model - stochastic linking and unlinking, (B) model 1 - inten- 
tional unlinking, (C) model 2 - preferential linking to, (D) model 3 - intentional 
unlinking and preferential linking to. The large frames present the general de- 
pendences on p and p; the inside frames show the change of magnetization if p 
is close to 0.5 and p is small. The lines without point markers are results for 
p = 0.001, 0.002, . . . , 0.01 subsequently. 




Fig. 2. Log plots of susceptibility of final configurations vs initial probability for 
a spin to be UP (p) for different values of short-cut edges (p) what is indicated 
by curve's labels: (A) model - stochastic linking and unlinking, (B) model 1 - 
intentional unlinking, (C) model 2 - preferential linking to, (D) model 3 - inten- 
tional unlinking and preferential linking to. The large frames present the general 
dependences on p and p; the inside frames show the change of susceptibility if 
p is close to 0.5 and p is small. The lines without point markers are results for 
p = 0.001, 0.002, . . . , 0.01 subsequently. 



- p E (0.48,0.52) in case of cellular automata with the network evolving ac- 
cording to the rules of models 0,1,2 

- p G (0.46,0.54) in case of cellular automata with the network evolving ac- 
cording to the rule of model 3. 

Because of the observed rapid changes in main characteristics one can say that 
the cellular automata solve the density classification task. Especially, model 
solves the density classification problem extremely efficiently: quickly and with 
the high certainty. 

Also time needed to reach the fixed point stabilization changes when p is 
increasing. If p > 0.1 and p « 0.5 then the time to stabilization is larger than 
in case of the system with not evolving network. However, the stabilization is 
reached in less than 40 steps. If p £ (0.4, 0.6) and p is small, p < 0.01, then the 
fixed point stabilization is not observed in less than 100 steps. The properties 
shown in figures collect features found as snap-shots of cellular automata states 
at the one hundredth step. This limitation does not introduce any important 
restriction. It is because that when the network does not evolve then, due to 
the many invariants of the spin-state rule, one observes the stabilizations on 
configurations which oscillate with some time period. This results in that the 
stabilization time is read as t = 100. Such limit configurations occur also if p 
is small, i.e., p < 0.01. When the rewiring process is stronger the oscillating 
patterns disappear in less than 100 steps. The stabilization different from the 
fixed points is not observed if p > 0.01. 

3.2 Degree distribution 

The vertex degree distributions are presented in Fig|3] The plots collect results 
recorded when the cellular automata systems arrive to the limit of either a ferro 
fixed point stabilization or a hundred step of time. Hence the total number of 
changes in the network depends on the time given to the evolution. However, 
if p > 0.01 then after a hundredth steps the average probability for an edge to 
change is 1. Thus, the difference between a network state after a hundred steps 
with p = 0.01 and after the first step with p = 1 consists in the synchronousness 
of events: with increasing p the edge evolution becomes more synchronous. 

All plots of the vertex degree distribution with p < 0.01 characterize the 
network configurations that are obtained at the one hundredth step. In case of 
model the fast convergence to the Gaussian distribution (the parabola shape in 
the log-plot) centered at k = 4 is observed. The significant distinction from the 
Gaussian shape in cases of other models is due to the preferences introduced. 
One observes vertices with the degree higher than 20 which in the Gaussian 
world do not exist. Because of the threshold value considered here, T = 8, the 
vertices with the degree larger than 8 are protected from unlinking and preferred 
in linking. Therefore the number of vertices with high degree is growing with 
increasing p. However, the presented distributions do not have wings depending 
polynomially on the vertex degree, namely, decaying as fc -7 for some 7 and k 
vertex degree . 




Fig. 3. Log plots of vertex degree in final configurations vs initial probability for 
a spin to be UP (p) for different values of short-cut edges (p) what is indicated 
by curve's labels: (A) model - stochastic linking and unlinking, (B) model 1 
- intentional unlinking, (C) model 2 - preferential linking to, (D) model 3 - 
intentional unlinking and preferential linking to. The gray lines are results for 
p = 0.001, 0.002, 0.009 subsequently; the dotted lines are results for p = 
0.2, 0.3, . . . , 0.9 subsequently. 



4 Concluding remarks 

We have considered properties of spin cellular automata initially spanned on 
the regular square lattice and then rewired each time step systematically. What 
we have observed is the final spin state. We have concentrated on the two con- 
figurations: all spins up and all spins down. It has appeared that only these 
two configurations are possible as the final configuration for spins. Which one of 
these two emerges is depended on p the density of 1 in the initial configuration. 
Due to the changes in the network the interval of uncertainty about the final 
configuration contracts sharply to the small interval around the critical value of 
p = 0.5. From the statistical physics point of view main investigations concern 
the ergodicity problem, i.e. the memory about the initial state properties in the 
evolving systems. The strong connectivity arisen from the small- world network 
enhances sharply the solution. The cellular automata considered are ergodic out- 
side the interval p G (0.48,0.52) what means that the interval of ergodicity is 
definitely larger than the interval of p in case of the ordinary spin cellular au- 
tomata. Therefore we claim that the systems considered can work as solvers of 
the density task. Moreover, since the analysed properties have been collected at 
time moments much smaller than the system size, the stated features seem to 
be universal in the sense that they are size independent. 

Our work is only preliminary. The dynamical system considered here needs 
further investigations. We have studied the zero-temperature Ising problem. Now 
the study of stochastic spin evolution should be undertaken. Moreover, the net- 
works considered here can bee included into the one thermodynamic ensemble, 
so-called canonical ensemble, because all networks have the same number of 
edges |16j . For the canonical ensemble there is a possibility to fix a temperature 
and associate the energy to every network configuration. Such approach states 
the new question: how this temperature and this energy relate to the ordinary 
thermodynamic characteristics of spin system. 
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